Jet Tomography of Quark Gluon Plasma 
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Recent experimental measurements of high pr hadron spectra and jet correlation at RHIC are 
analyzed within a parton model which incoporates initial jet production and final propagation 
in heavy-ion collisions. The suppresion of single hadron spectra, back-to-back correlation, their 
centrality dependence and azimuthal anisotropy point to a dense matter with an initial parton 
density about 30 times of that in a cold heavy nucleus. 



I. INTRODUCTION 

In high-energy heavy-ion collisions, a dense medium 
of quarks and gluons is expected to be produced and 
possibly a quark-gluon plasma is formed. One impor- 
tant step in establishing evidence of QGP formation is to 
charaterize the properties of the dense medium produced, 
for example, the parton and energy density and color 
confinement, among many other charateristics. Conven- 
tionally, one can study the properties of a medium via 
scattering experiments with particle beams. In deeply 
inelastic scattering (DIS) experiments, for example, lep- 
tons scatter off the nucleon medium via photon exchange 
with quarks. The response function or the correlation 
function of the electromagnetic currents, 
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is a direct measurement of the quark distributions in a 
nucleon or nucleus. Such experiments have provided un- 
precedented information about partonic structure of nu- 
cleons and nucleus and confirmed the prediction of QCD 
evolution [1]. 

For dynamic systems such as that produced in heavy- 
ion collisions, one can no longer use the technique of scat- 
tering with a beam of particles because of the transient 
nature of the matter. The lifetime of the system is very 
short, on the order of a few fm/c. The initial spatial size 
is only the size of the heaviest nuclei, about 6 fm in di- 
ameter in the transverse dimension. The system expands 
very rapidly both in the longitudinal and transverse di- 
rection. These characteristics make it impossible to use 
external probes to study the properties of the produced 
dense matter in high-energy heavy- ion collisions. For- 
tunately, one can prove that the thermal average of the 
above correlation function gives the photon emission rate 
from the evolving system. The emission rate depends 
mainly on the local temperature or the parton density 
while the total yield also depends the whole evolution 
history of the system. Therefore, a strongly interacting 
system can reveal its properties and dynamics through 
photon and dilepton emission. One can further study 
the resonance properties of the emitted virtual photons 
and their medium modification. The screening of strong 
interaction in a color deconfincd medium should lead to 



dissociation of the binding states and thus the quarko- 
nia suppression [2]. Such color screening is a result of 
strong interaction between quarks and gluons at high 
density and temperature. The same interaction will also 
cause attenuation of fast and energetic partons propagat- 
ing through the medium. Such an effect is the underly- 
ing physics of the jet quenching [3] phenomenon and jet 
tomography technique for studying properties of dense 
matter in high-energy heavy-ion collisions. 

Jet quenching as a probe of the dense matter in heavy- 
ion collisions, takes advantage of the hard processes of 
jet production in high-energy heavy-ion collisions. Sim- 
ilar to the technology of computed tomography (CT), 
study of these energetic particles, their initial production 
and interaction with the dense medium, can yield criti- 
cal information about the properties of the matter that 
is otherwise difficult to access through soft hadrons from 
the hadronization of the bulk medium. Though relatively 
rare with small cross sections, the jet production rate can 
be calculated perturbative in QCD and agrees well with 
experimental measurements in high-energy pp(p) colli- 
sions. A critical component of the jet tomography is then 
to understand the jet attenuation through dense matter 
as it propagates through the medium. 



II. MODIFIED FRAGMENTATION FUNCTION 

A direct manifest of jet quenching is the modification 
of the fragmentation function of the produced parton, 
D a ^h(z, /j?) which can be measured directly. This mod- 
ification can be directly translated into the energy loss of 
the leading parton. 

To demonstrate medium modified fragmentation func- 
tion and parton energy loss, one can study deeply inelas- 
tic scattering (DIS) eA [4-6]. Here, we consider the semi- 
inclusive processes, e(Li) + A(p) — > e(L 2 ) + h(£h) + X, 
where L± and L 2 are the four-momenta of the incoming 
and the outgoing leptons, and £h is the observed hadron 
momentum. The differential cross section for the semi- 
inclusive process can be expressed as 



a 



EM 



1 



d 3 L 2 d 3 £ h 2tts Q 



-L fJi „Ee h - 



~dHZ 



(2) 



where p = [p + ,0, 0j_] is the momentum per nucleon in 
the nucleus, q = L 2 — L\ = [— Q 2 /2q~, q~ , 0j_] the mo- 
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mentum transfer, s = (p + Li) 2 and oem is the electro- 
magnetic (EM) coupling constant. L M „ is the leptonic 
tensor while W^ v is the semi-inclusive hadronic tensor. 

In the parton model with the collinear factorization ap- 
proximation, the leading-twist contribution to the semi- 
inclusive cross section can be factorized into a product of 
parton distributions, parton fragmentation functions and 
the hard partonic cross section. Including all leading log 
radiative corrections, the lowest order contribution from 
a single hard 7* + q scattering can be written as 



dzh 



where Hffi (x,p, q) is the hard part of the process in lead- 
ing order, the momentum fraction carried by the hadron 
is defined as Zh = ^Z/<1~ and xb = Q 2 /2p + q~ is the 
Bjorken variable. \x\ and /x 2 are the factorization scales 
for the initial quark distributions /^(x,^ 2 ) in a nucleus 
and the fragmentation functions D q ^h(zh, M 2 ) > respec- 
tively. 

In a nuclear medium, the propagating quark in DIS will 
experience additional scatterings with other partons from 
the nucleus. The rescatterings may induce additional 
gluon radiation and cause the leading quark to lose en- 
ergy. Such induced gluon radiations will effectively give 
rise to additional terms in the evolution equation leading 
to the modification of the fragmentation functions in a 
medium. These are the so-called higher-twist corrections 
since they involve higher-twist parton matrix elements 
and are power-suppressed. We will consider those con- 
tributions that involve two-parton correlations from two 
different nucleons inside the nucleus. 

One can apply the generalized factorization to these 
multiple scattering processes [7]. In this approximation, 
the double scattering contribution to radiative correction 
can be calculated and the effective modified fragmenta- 
tion function is 
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D q ^ h {z h ,n ) = D q ^ h {z h ,n ) + -jr^ J — 

[A<y q ^ qg (z,x,x L ,e^) D q ^ h (z h /z) 
+ A- /q ^ gq (z,x,x L , l%)D g _ h {z h /z)] , (4) 

where D q ^h{zh, M 2 ) and D g ^h{zh, M 2 ) arc the leading- 
twist fragmentation functions. The modified splitting 
functions are given as 
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Here, the fractional momentum is defined as ii = 
l\j2p+q~ z(l - z) and x = x B = Q 2 /2p + q~ is the 



Bjorken variable. The twist-four parton matrix elements 
of the nucleus, 
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has a dipolc-like structure which is a result of 
Landau-Pomerachuck-Migdal (LPM) interference in 
gluon brcmsstrahlung. Here the intrinsic transverse mo- 
mentum is neglected. In the limit of collinear radiation 
(xl — ► 0) or when the formation time of the gluon radia- 
tion, Tf = 1/xlp + , is much larger than the nuclear size, 
the destructive interference leads to the LPM interference 
effect. 

Using the factorization approximation [4,5,7,8], we can 
relate the twist-four parton matrix elements of the nu- 
cleus to the twist-two parton distributions of nucleons 
and the nucleus, 
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where C = 2CxTfg{xr) is considered a constant. One 
can identify 1/xlp + = 2q~ z(\ — z)jl\ as the formation 
time of the emitted gluons. When it becomes comparable 
or larger than the nuclear size, the above matrix element 
vanishes, demonstrating a typical LPM interference ef- 
fect. 

Since the LPM interference suppresses gluon radia- 
tion whose formation time (77 ~ Q 2 /&rP + ) is larger 
than the nuclear size MRa/p + in our chosen frame, £ T 
should then have a minimum value of £ T ~ Q 2 /MRa ~ 
Q 2 /A 1 / 3 . Here M is the nucleon mass. Therefore, 
the leading higher-twist contribution is proportional to 
asRA/Bp <~ a s R 2 A /Q 2 due to double scattering and de- 
pends quadratically on the nuclear size Ra ■ 

With the assumption of the factorized form of the 
twist-4 nuclear parton matrices, there is only one free 
parameter C(Q 2 ) which represents quark-gluon correla- 
tion strength inside nuclei. Once it is fixed, one can pre- 
dict the z, energy and nuclear dependence of the medium 
modification of the fragmentation function. Shown in 
Fig. 1 are the calculated nuclear modification factor of 
the fragmentation functions for 14 N and 84 Kr targets 
as compared to the recent HERMES data [9]. The pre- 
dicted shape of the z-dependence agrees well with the 
experimental data. A remarkable feature of the predic- 
tion is the quadratic A 2 / 3 nuclear size dependence, which 
is verified for the first time by an experiment. By fitting 
the overall suppression for one nuclear target, we obtain 
the only parameter in our calculation, C(Q 2 ) = 0.0060 
GeV 2 with a s (Q 2 ) = 0.33 at Q 2 w 3 GeV 2 . The pre- 
dicted ^-dependence also agrees with the experimental 
data [22]. 
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FIG. 1. Predicted nuclear modification of jet fragmenta- 
tion function is compared to the HERMES data [9] on ratios 
of hadron distributions between A and d targets in DIS. 

We can quantify the modification of the fragmentation 
by the quark energy loss which is defined as the momen- 
tum fraction carried by the radiated gluon, 
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where x, 



p 2 /2p + q~ z(l — z) = xb/z(1 — z) if we choose 
the factorization scale as p 2 = Q 2 . When x A <C xb <C 1 
we can estimate the leading quark energy loss roughly as 
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Since x A — 1/MR A , the energy loss (Az g ) thus depends 
quadratically on the nuclear size. 

In the rest frame of the nucleus, p + = to at, q~ = 
v, and x B = Q 2 /2p + q~ = Q 2 j2m^v . One can get 
the averaged total energy loss as AE = v(Az g ) w 
C(Q 2 )a 2 (Q 2 )m N R\(C A /N c )3\n(l/2x B ). With the de- 
termined value of C, (x B ) w 0.124 in the HERMES ex- 
periment [9] and the average distance (La) = R A \j2j-K 
for the assumed Gaussian nuclear distribution, one gets 
the quark energy loss dE/dL w 0.5 GeV/fm inside a Au 
nucleus. 



III. JET QUENCHING IN HOT MEDIUM AT 
RHIC 

To extend our study of modified fragmentation func- 
tions to jets in heavy-ion collisions, we can assume (fcf.) w 



/i 2 (the Debye screening mass) and a gluon density pro- 
file p(y) — (tq/t)9(R a — y)po for a 1-dimensional ex- 
panding system. Since the initial jet production rate 
is independent of the final gluon density which can be 
related to the parton-gluon scattering cross section [10] 
[a s XTG(xT) ~ /i 2 (Tg], one has then 
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where r/ = 2Ez{\ — z)/£?n is the gluon formation time. 
One can recover the form of energy loss in a thin plasma 
obtained in the opacity expansion approach [11], 
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Keeping only the dominant contribution and assuming 
C a 2na 2 ./p 2 (C a =l for qg and 9/4 for gg scattcr- 



cr, 



ing), one obtains the averaged energy loss, 



dE _ TTC a C A al 
{ dL } ~ R A 



pRa 

I drp{ 

J Tn 



r)(r 



t ) In — -. 

Tp, Z 



(13) 



Neglecting the logarithmic dependence on r, the aver- 
aged energy loss in a 1-dimcnsional expanding system 
can be expressed as (4f )id ~ (dE^/ dL)(2ro/ Ra), where 
dE /dL oc poR A is the energy loss in a static medium 
with the same gluon density po as in a 1-d expanding 
system at time To. Because of the expansion, the aver- 
aged energy loss (dE/dL) id is suppressed as compared 
to the static case and does not depend linearly on the 
system size. 

In order to calculate the effects of parton energy loss 
on the attenuation pattern of high pt partons in nuclear 
collisions, we use a simpler effective modified fragmenta- 
tion function [12,13], 
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where z' c , z' g are the rescaled momentum fractions. This 
effective model is found to reproduce the pQCD result 
from Eq.(4) very well, but only when Az = AE C /E is 
set to be Az Q.6(z g ). Therefore the actual averaged 
parton energy loss should be AE/E — 1.6Az with Az 
extracted from the effective model. The factor 1.6 is 
mainly caused by the unitarity correction effect in the 
pQCD calculation. 

Since gluons are bosons, there should also be stimu- 
lated gluon emission and absorption by the propagating 
parton because of the presence of thermal gluons in the 
hot medium. Such detailed balance is crucial for parton 
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thermalization and should also be important for calculat- 
ing the energy loss of an energetic parton in a hot medium 
[14]. Taking into account such detailed balance in gluon 
emission, one can then get the asymptotic behavior of the 
effective energy loss in the opacity expansion framework 
[14], 
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where the first term is from the induced bremsstralung 
and the second term is due to gluon absorption in detailed 
balance which effectively reduce the total parton energy 
loss in the medium. Numerical calculations show that 
the effect of the gluon absorption is small and can be ne- 
glected for partons with very high energy. However, the 
thermal absorption reduces the effective parton energy 
loss by about 30-10% for intermediate values of parton 
energy. This will increase the energy dependence of the 
effective parton energy loss in the intermediate energy 
region. One can parameterize such energy dependence 
as, 



, dE . . . 



1.6) 1 - 2 /(7.5 + £//i), 



(16) 



The threshold is the consequence of gluon absorption that 
competes with radiation that effectively shuts off the en- 
ergy loss. The parameter p is set to be 1 GeV in the 
calculation. 

To calculate the modified high pt spectra in A + A 
collisions, we use a LO pQCD model [15,16], 
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with medium modified fragmentation funcitons D' h ^ c 
given by Eq. (14) and the fragmentation functions in free 
space D^ c (z c ,Q 2 ) are given by the BBK parameteriza- 
tion [17]. Here, z c = pt/ptc, y = Ucj &{ab cd) are ele- 
mentary parton scattering cross sections and t A {b) is the 
nuclear thickness function normalized to J d 2 bt A (b) = A. 
We will use a hard-sphere model of nuclear distribution 
in this paper. The K « 1.5 — 2 factor is used to account 
for higher order pQCD corrections. 

The parton distributions per nucleon / a M(x a , Q 2 , r) 
inside the nucleus are assumed to be factorizable into 
the parton distributions in a free nucleon given by 
the MRSD— ' parameterization [18] and the impact- 
parameter dependent nuclear modification factor which 
will given by the new HIJING parameterization [19]. The 



initial transverse momentum distribution g A (kx, Q 2 , b) is 
assumed to have a Gaussian form with a width that in- 
cludes both an intrinsic part in a nucleon and nuclear 
broadening. This model has been fitted to the nuclear 
modification of the px spectra in p + A collisions at up 
to the Fermilab energy -Js = 40 GeV. Shown in Fig. 2 
are the first prediction made in 1998 [15] of the Cronin 
effect at RHIC for p + Au collisions at y/s — 200 GeV as 
compared to the recent RHIC data. As one can see, the 
initial multiple scattering in nuclei can give some moder- 
ate Cronin enhancement of the high px spectra. There- 
fore, any suppression of the high px spectra in Au + Au 
collisions has to be caused by jet quenching. 
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FIG. 2. The first predictions [15] of the Cronin effect in 
p + Au collisions at yfs — 200 GeV are compared to the re- 
cent RHIC data from PHENIX [20] and STAR [21]. 

We assume a 1-dimensional expanding medium with a 
gluon density p g (r, r) that is proportional to the trans- 
verse profile of participant nucleons. According to Eq. 13, 
we will calculate impact-parameter dependence of the en- 
ergy loss as 
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(18) 



where AL(&, r, <p) is the distance a jet, produced at r, 
has to travel along ft at an azimuthal angle <f> relative to 
the reaction plane in a collision with impact-parameter 
b. Here, po is the averaged initial gluon density at tq 
in a central collision and (dE/dL)id is the average par- 
ton energy loss over a distance R A in a 1-d expand- 
ing medium with an initial uniform gluon density p . 
The corresponding energy loss in a static medium with 
a uniform gluon density po over a distance R A is [22] 
dE /dL — (R A /2ro)(dE /dL)id. We will use the parame- 
terization in Eq. (16) for the effective energy dependence 
of the parton quark energy loss. 
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Shown in Fig. 3 are the calculated nuclear modifica- 
tion factors Rab(pt) = da\ B J (A/binary) d(Tp v for hadron 
spectra (\y\ < 0.5) in Au + Au collisions at y/s = 200 
GeV, as compared to experimental data [24,23]. Here, 

(Aginary) = J d 2 bd 2 rt A (r)t A (\b - r|). To fit the observed 
7T° suppression (solid lines) in the most central collisions, 
we have used fi — 1.5 GeV, e = 1.07 GeV/fm and 
Aq = l/(<jpo) = 0.3 fm. The hatched area (also in other 
figures in this paper) indicates a variation of eo = ±0.3 
GeV/fm. The hatched boxes around Rab = 1 repre- 
sent experimental errors in overall normalization. Nu- 
clear kx broadening and parton shadowing together give 
a slight enhancement of hadron spectra at intermediate 
Pt = 2 — 4 GeV/c without parton energy loss. 

The flat px dependence of the ir° suppression is a con- 
sequence of the strong energy dependence of the parton 
energy loss. The slight rise of Rab at px < 4 GeV/c 
in the calculation is due to the detailed balance effect 
in the effective parton energy loss. In this region, one 
expects the fragmentation picture to gradually lose its 
validity and is taken over by other non-perturbative ef- 
fects, especially for kaons and baryons. As a consequence, 
the (K + p)/tt ratio in central Au + Au collisions is sig- 
nificantly larger than in peripheral Au + Au or p + p 
collisions. To take into account this effect, we add a 
nuclear dependent (proportional to (Aginary)) soft com- 
ponent to kaon and baryon fragmentation functions so 
that (K + p)/n ~ 2 at px ~ 3 GeV/c in the most cen- 
tral Au + Au collisions and approaches its p + p value 
at pt > 5 GeV/c. The resultant suppression for total 
charged hadrons (dot-dashed) and the centrality depen- 
dence agree well with the STAR data. One can directly 
relate and it suppression via the (K + p) /tt ratio: 



nh ± 
n AA 



R AA [1 + {K + P )/tt]aa/[1 +{K+ p)/n] pp . It is 



clear from the data that (K + p) /it becomes the same 
for Au + Au and p + p collisions at px > 5 GeV/c. To 
demonstrate the sensitivity to the parameterized parton 
energy loss in the intermediate px region, we also show 
Raa m 0-5% centrality (dashed line) for a = 2.0 GeV 
and eo = 2.04 GeV/fm without the soft component. 

In the same LO pQCD parton model, one can also 
calculate di-hadron spectra, 
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(19) 



for two back-to-back hadrons from independent fragmen- 
tation of the back-to-back jets. Let us assume hadron hi 
is a triggered hadron with pxi = Px lg - One can de- 
fine a hadron-triggered fragmentation function (FF) as 



the back-to-back correlation with respect to the triggered 
hadron: 



D h ^{z Tl 4>J^)=p t ^ 



da^ 2 /d 2 p^ ls dp T d(l) 
da h A \/d 2 p^ 



2 trig 



(20) 



similarly to the direct-photon triggered FF [12,13] in 7- 
jet events. Here, zt — Pt/Pt 1S and integration over 
12/1,2 1 < Ay is implied. In a simple parton model, the 
two jets should be exactly back-to-back. The initial par- 
ton transverse momentum distribution in our model will 
give rise to a Gaussian-like angular distribution. In ad- 
dition, we also take into account transverse momentum 
smearing within a jet using a Gaussian distribution with 
a width of (k±) = 0.6 GeV/c. Hadrons from the soft 
component are assumed to be uncorrelated. 

Shown in Fig. 4 are the calculated back-to-back cor- 
relations for charged hadrons in Au + Au collisions as 
compared to the STAR data [25]. The same energy 
loss that is used to calculate single hadron suppres- 
sion and azimuthal anisotropy can also describe well the 
observed away-side hadron suppression and its central- 
ity dependence. In the data, a background B(j>t)[1 + 
2v 2 (pt) cos(2A0)] from uncorrelated hadrons and az- 
imuthal anisotropy has been subtracted. The value of 
v%{pr) is measured independently while B(j>t) is deter- 
mined by fitting the observed correlation in the region 
0.75 < |0| < 2.24 rad [25]. 
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FIG. 3. Hadron suppression factors in Au + Au collisions 
as compared to data from STAR [24] and PHENIX [23]. See 
text for a detailed explanation. 
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FIG. 4. Back-to-back correlations for charged hadrons with 
Py ig > p T > 2 GeV/c, p* is = 4-6 GeV/c and |y| < 0.7 in 
Au + Au (lower curves) and p + p (upper curves) collisions as 
compared to the STAR [25] data. 

With both the single spectra and dihadron spectra, 
the extracted average energy loss in this model calcu- 
lation for a 10 GeV quark in the expanding medium is 
{dE/dL) ld ps 0.85 ± 0.24 GeV/fm, which is equivalent 
to dEo/dL ps 13.8 ± 3.9 GeV/fm in a static and uniform 
medium over a distance Ra = 6.5 fm. This value is about 
a factor of 2 larger than a previous estimate [22] because 
of the variation of gluon density along the propagation 
path and the more precise RHIC data considered . 



IV. SUMMARY 

In summary, with the recent measurements of high px 
hadron spectra and jet correlations in d+Au collisions at 
RHIC, the observed jet quenching in Au + Au collisions 
has been established as a consequence of final state in- 
teraction between jets and the produced dense medium. 
The collective body of data: suppression of high px spec- 
tra and back-to-back jet correlation, high pt anisotropy 
and centrality dependences of the observables, indicate 
that the cause of the jet quenching is due to parton en- 
ergy loss rather than final state absorption of hadrons. 

A simultaneous phcnomcnological study of the sup- 
pression of hadron spectra and back-to-back correla- 
tions, and high-pT azimuthal anisotropy in high-energy 
heavy-ion collisions within a single LO pQCD parton 
model incorporating current theoretical understanding 
of parton energy loss can describe the experimental 
data of Au + Au collisions very well. With HIJING 
(EKS) parton shadowing, the extracted average energy 



loss for a 10 GeV quark in the expanding medium is 
(dE/dL) ld ps 0.85(0.99) ± 0.24 GeV/fm, which is equiv- 
alent to dE /dL ps 13.8(16.1) ± 3.9 GeV/fm in a static 
and uniform medium over a distance Ra = 6.5 fm. This 
value is about a factor of 2 larger than a previous esti- 
mate [22] because of the variation of gluon density along 
the propagation path and the more precise RHIC data 
considered here . 
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